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The transient network literature up to now has considered that the connection probabihty of a 

, free strand does not depend on the strand extension, in contrast with the disconnection probabihty. 

' We argue that, on thermodynamic grounds, both probabihties must have the same dependence on 

«' I the strand extension. We propose a model for the pinning/snapping of elastic strands reversibly 

. binding to sticky sites, giving explicit expressions for the connection and disconnection rates. 
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INTRODUCTION 



This paper addresses a seemingly paradoxical result uncovered when assessing the thermodynamic consistency 
of Yamamoto's transient network model (Yamamoto 1956,1957,1958, Tanaka and Edwards 1992a-d, Palierne 2000): 
Consider an elastic strand, i.e. a chainlike molecule or part of molecule, with both ends pinned to fixed sticky sites by 
O , reversible bounds. It is fairly obvious that the probability per unit time that a pinned end snaps back and becomes 
^ . free is an increasing function of the strand extension, hence of the separation of the sticky sites. It is less obvious, 
'"^ ' however, that the probability per unit time that a free strand becomes pinned when encountering a sticky site should 
, increase in the same proportion with the strand extension. This fact stems from the detailed balance condition for 
^ ' the connection/disconnection reaction formulated in Palierne (2000), the violation of this condition being at the root 
of the thermodynamic inconsistencies in Yamamoto's model and in its offspring (Palierne 2000). The argument is as 
follows. 

G [ We take one end of a strand fixed at the origin of space, and consider the other end. Its position is h, and we 
' 1 consider its pinning to a sticky site, a process hereafter called the connection of a free strand, as well as the reverse 
^ ' process called disconnection, whereby the pinned end snaps back and the connected strand becomes free. Connection 
Q , and disconnection are instantaneous and take place without change of the connector h. (hereafter, notations close 
Q ' to those of Palierne 2000 will be used, except for introducing the one-strand distribution ipa{ii) — N^^Tah, where 
a = F or C denotes the connection state). In the ideal systems we consider, the strand chemical potential /iQ(h) is 
^ ' the sum of the distribution entropy kBTlnipaO^) and a part CQ(h) depending on the number and the energy of the 
^ strand conformations compatible with the h value of the connector. The chemical potential at equilibrium is a— and 
h— independent, /^^(h) = fi'^'', and the equilibrium distribution reads 



a=F,C 



V'o(h)=ife-^°('')/'=-^ , K-'^ V e-'^-^''^^'-^ d^h (1) 
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' where K ensures the normalisation / {ipc{i>-) + '0F(h))(i'^h = 1. The h-dependence of Caih) only involves the configu- 

■>«.^ ration of the strand, irrespective of its connection state, therefore the connection potential, defined as the difference 

AC = Cc{h)-CF{h), (2) 

I ' can not depend on h, and V'c(h) must be proportional to iljp{h.): ipcO^) = e~^''/'^^-^V'i="(h) For instance, freely-jointed 

^5 . strands comprising n segments of length a have C-F(h) = Cc(h) — = {SksT / 2na^)h^ and both the free and the 

O connected strands are Gaussian: V'c(h) oc ippih) cx e~^^ /^na (Paherne 2000). 

O I Consider now the kinetics of the connection/disconnection reaction. Let /3c (h) denote the connection rate, i.e. the 

^ ' probability per unit time that a free strand connects, and let accordingly (3p{h) denote the disconnection rate. At equi- 

• i-H , librium, the distributions V'-F(h) and ipcO^) of free and connected strands are such that the connection/disconnection 
reaction rate 

. ; ^(h) = /3c(h)^F(h) - PFmdh) (3) 



vanishes, resulting in the detailed balance condition 

/3c (h) _ jjcjh) 
l3F{h) VF(h) 



g-AC/fcsT independent ofh (4) 



The inescapable consequence of this relation is, if we admit that the disconnection rate /3p{h) increases with h, then 
the connection rate /3c (h) must increase in the same proportion. In the following section, both /3_F(h) and /3c (h) will 
be derived from a microscopic model, and relation ^ will be checked. 
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MICROSCOPIC MODEL 



We are interested in the interaction between a sticky site and one strand's end we call the sticker. The other end 
will therefore be kept fixed at the origin of space. Let {c. h} denote the set of degrees of freedom, h being the position 
of the non-fixed end and c denoting the other degrees of freedom, c and h are assumed to vary independently and to 
fluctuate according to an equilibrium Gibbs distribution at temperature T. The energy of configuration {c, h} is the 
sum of the strand internal energy e(c, h) and the interaction energy 

es{h)=we{a-\h.-s\) (5) 

describing a spherical well centered on s, of depth w and radius a, attracting the sticker and no other part of the 
strand (w < for an attracting well). The link between this microscopic model and macroscopic thermodynamics is 
provided by the quantity 

C(h) = -/cBTln^e-^'^'B-"'' = Co(h) + e4h) (6) 



where Co(h) = -feTln^ e ^ (7) 

C 

where designates teh statistical sum over all configurations keeping a given h. The conector force 

f(h) = ^ (8) 

is the statistical average of the microscopic force 9e(c, h)/9h over the Gibbs distribution of c. The radius a will be 
considered small with respect to the length scale of e, say a <C \Je/V^e, but the product af(h) need not be smaller 
than UbT. Expanding Co(h) about s according to Co(h) = Co(s) + f(s) • (h — s), we write down for further use the 
following integrals 

Io= e (fh (9) 

f <:o(h) co(3) 

/(s) = / e "bt d^h = Awa^e "^(w"^ coshw - u"^ sinhw) (10) 

^|h-s|<o 

r Co(s+a) Cn(s) , 

J(s) = / e '=bt d^a = 47ra^e "bt sinhu (11) 



f Co(h) + es(h) 

K{s)= e ^ d^h = Jo + (e"'^ - l)/(s) (12) 

where u = a\i{s)\/kBT. /o and K{s) are integrals over the whole h space, /(s) over the sphere |h — s| < a, and J(s) 
over the surface of the same sphere. The equilibrium Gibbs distribution of h reads 



Co(h)+.a(h) 
k„T 



K{s) 

The probability p of finding the sticker within the well can be written 



(13) 



We now introduce the kinetic properties of our model. The sticker is assumed to diffuse in a Brownian fashion both 
within and outside the well. The surface |h — s| = a is assumed to be a potential barrier such that Aout) the probability 
per unit time and per unit area that the sticker crosses the barrier from the inside to the outside at point s + a reads 
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Aout(s + a) = tout(/'in(s + a) = tout lim 0(s + (1 + 7/)a), (15) 

tout being the outwards transmission factor, and the probabihty of the reverse crossing can be written accordingly 

Ain(s + a) = ii„0out(s + a) = iin Urn (j){s + {1 + T])a), (16) 



tout and tin are intrinsic properties of the barrier, independent of the position on the sphere. According to eq. 13, the 
hmits (t>in and 0out are such that 

pt±^ = (17) 
0out(s + a) 

independent of the position on the sphere. The principle of detailed balance states that 

Ain — Aout (18) 
at equilibrium, so the transmission factors are linked to the well depth by the customary relation 

— = e"^ (19) 

^out 

The escape probability per unit time of a sticker in a well reads then 

1 /" I „^J2 _ toutJ{s) _ tout u 



ms) = - / Aout(s + a)d^a= = — TT^ (20) 

Pj\a\=a ^(s) a C{u) 

where u — a|f (s)|/A;Br, and C{u) = cothw — u^^ is the Langevin function. One has u/C{u) = 3 + + ... for u ^ 
and u/C{u) — u + 1 + o(l) for u ^ cx) ; and u/ C{u) = (1 + 2(1 + /4)^) over the whole range of u to within 3% 
precision. 



DISTRIBUTION OF STRANDS AND WELLS 

We model the network as a system of non-interacting strands in a statistically homogeneous distribution of attracting 
wells. The wells number density n is such that many wells enter the volume spanned in the sticker's Brownian 
motion. The results of the preceding section need but little modification: using the identity J I{s)d^s — ^Tra^Io, the 
normalisation integral can be rewritten 

K = Io + (e^ - l)n J /(s)rf^s = (1 - x + xe'^)Io (21) 

where x — ^irna^ is the volume fraction occupied by the wells. The distribution for a given strand becomes 

ie ''B^ if h outside a well (probability 1 — x) 

f _^ (22) 
^gi-BTg kgT if h inside a well (probability x) 

The ensemble of strands is thus described by the connected strand distribution 

X Co(h) 

V'c(h) = — e^B^e "B-r (23) 
K 



and the free strand distribution 



1 — X _ioM 

V'F(h) = -^^e "bt (24) 
K 



satisfying the one-strand normalisation of eq. (0), J (V'c(h) + ■(/'i='(h))d'^h — 1. The ratio 

■0c (h) X 



I (x.\ 1 ^'^^ (25) 

TpFK'^) l — X 
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takes the form (Q) upon introducing the connection chemical potential 

AC^w + fcaTlni— ^ (26) 

X 

as the sum of the energy w and the distribution entropy fceTln -!— ^ accounting for the different volume fraction offered 
to the connected and to the free stickers. The present model thus satisfies the equilibrium properties announced in 
the introduction. 

The disconnection rate /?_F(h) is the escape probability per unit time for a connected strand, i.e. a strand whose 
sticker lies in a well centered on s such that |h — s| < a. At lowest order in a, we simply have Vl{s) — ri(h), and 

/3f^(h) = il7(h) (27) 

X 

where the fraction 1/x compensates for the normalisation factor x appearing in cq. (2^), so that l3F{^)4'ci)^) = 



ri(h)0(h in a well). The conection rate /3c(h) relates the probability per unit time of fa 
strand distribution, so one has similarly 



ling into the well to the free 



/3c(h) = 0(h) (28) 

1 — x 

It must be noted that our I3pi)i) has a much weaker dependence in h than the exponential form proposed by Tanaka 
and Edwards (1992b), whereas their /3c(h) is a constant, resulting in a connected strand distribution (their eq. 2.23) 
incompatible with the Gibbs distribution. 

Our /3i?(h) and /3c (h) are found to satisfy the detailed balance (^: 

Ml) - f29) 
In fact, relation(|3) expresses the detailed balance at the global level of connection/disconnection reaction, whereas 



relation (29) is a consequence of the detailed balance at the local level of crossing the barrier at the well border. 
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